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Abstract 

A methodology  is  presented  for  analyzing  periodic  autoregressions 
which  is  also  applicable  when  inferring  the  second  order  properties  of 

periodically  correlated  processes.  In  addition,  capitalizing  on  the  duality 

between  periodic  and  multiple  autoregressions,  a method  is  set  forth  for 

analyzing  the  latter,  which  overcomes  the  usual  requirements  of  a large 

number  of  both  parameters  and  computer  storage  locations. 
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1.  Introduction 

Given  a time  series  [Y(t)  ; t ■ 0,±1,...]  whose  second  order  moments 
exist,  define  its  mean  function,  m(t)  ° EY( t)  , and  its  covariance  kernel, 
R(s,t)  ■ e|y(s)  - m(s)J  ^Y(t)  - m(t5j  • A class  of  nonstationary  processes 
which  readily  lends  itself  to  analysis,  and  is  of  practical  importance 
(see  [9],  [5],  [18])  is  that  of  periodically  correlated  processes  (see  [2]). 

Definition.  The  process  Y(*)  is  said  to  be  periodically  correlated  of 
period  d , if  for  some  positive  integer  d and  for  all  Integers  s,t  , 

m(t)  ■ m(t  + d)  , R(s,t)  ■ R(s  + d,t  + d) 

Since  we  propose  dealing  with  the  second  order  properties  of  the 
process,  without  loss  of  generality  take  m(t)  ■ 0 . 

Periodically  correlated  processes  are  not  only  of  interest  in  their 
own  right,  but,  because  of  their  duality  with  multivariate  covariance 
stationary  time  series,  they  also  provide  insight  into  and  modeling  facility 
for  these  series.  This  claim  is  based,  in  the  main,  upon  the  following  con- 
struction: define  the  Jth  component  of  the  d-dimensional  vector  X(t)  by 

(see  [2]), 

(1.1)  Xj(t)  - YQ+d(t-ljJ 
and  the  covariance  kernel  of  X(»)  by 

Rjk(s,t)  - EXj(s)  Xfc(t)  , 1 s ],k  s d , s,t-0,±l,... 


f 

A 
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By  noting  that 

(1.2)  RJk(s,t)  - R(j+ds,k+dt)  , 
we  nave, 

Theorem  (Gladyshev  [2]).  The  process  Y(*)  is  periodically  correlated  of 
period  d if,  and  only  if,  the  process  X(*)  is  covariance  stationary. 

(V 

This  theorem  shows  that,  although  not  covariance  stationary,  the 
process  Y(*)  , nonetheless,  does  not  deviate  too  much  from  stationarity ; 
and,  in  fact,  its  second  order  properties  may  be  deduced  from  those  of 
X(*)  . We  can  further  capitalize  on  this  duality  but  in  the  other  direction; 
given  a covariance  stationary  process  X(*)  define  the  associated  periodi- 
cally correlated  process,  Y(*)  , by  (1.1).  This  is  especially  useful  when 
X(*)  is  a multiple  autoregression.  For  then  by  investigating  Y(«)  , we 
actually  effect  an  easily  calculated  orthogonal  decomposition  of  the  process 
£(*)  > thus  achieving  the  usual  gain  in  simplicity  and  power  associated  with 
orthogonalization.  This  point  is  amplified  in  the  remainder  of  the  paper. 
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2.  Autoregressions 

Given  a sample  X(l) , . . . ,X(T)  from  a d-dimensional  zero  mean  covar- 
iance  stationary  time  series  with  spectral  density  matrix  f(w)  , 

-rr  s U)  s;  rr  , the  problem  is  to  estimate  f(»)  . Under  some  very  mild 
restrictions  on  f(*)  , [7],  it  can  be  written  as  the  spectral  density  of 
an  infinite  order  autoregression.  Parzen  [13]  thus  proposes  treating  the 
sample  as  one  from  a £th  order  autoregression, 

P 

(2.1)  X(t)  + E A(j)  X(t  - j)  - e(t) 

j-l 

where  cov(e(t))  ■ E , and  by  choosing  p large  enough,  f(*)  can  be 
arbitrarily  closely  approximated.  Parzen  [15]  gives  a method  for  choosing 
p which  accomodates  both  the  numerical  approximation,  which  argues  for  a 
large  p , and  the  requisite  statistical  estimation,  which  argues  for  a 
small  p . The  autoregressive  approximants  have  been  used  in  practice 
with  success,  (see  [4],  [14]  and  references  therein,  also  in  geophysics 
[16]  where  a closely  analogous  method,  the  method  of  maximal  entropy,  is 
used).  For  these  reasons  it  seems  Important  to  study  autoregressive 
processes. 

The  periodically  correlated  analogue  of  an  autoregression  is  given 
by  (see  [5]), 

Definition.  A process  Y(*)  is  said  to  be  a periodic  autoregression  of 
period  d and  order  (p^,...,p(j)  if  for  all  integer  t , 

pt 

(2.2)  Y(t)  + E a (J)  Y(t  - j)  - fi(t) 

J-l 


‘•5* 
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2 2 

where  the  e(*)  are  uncorrelated  with  mean  zero  and  Ee  (t)  - at  > and, 
pt  " Pt-Kl  » CTt  " at+d  » and*  at<J)  " at-hi(J)  * i " * 


Theorem  1.  If  X(*)  and  Y(*)  are  associated  by  (1.1),  then,  X(*)  is 

an  autoregression  of  order  p with  positive  definite  2 if,  and  only  if, 

Y(*)  is  a periodic  autoregression  of  period  d and  order  (p^,...,Pd) 

2 2 

with  positive  n1>.*«»Od  > and, 

p - [max  (p  - j)/d]  + 1 , 

3 J 

where,  for  integral  j , [xl  - j for  j £ x < j + 1 . 


Proof.  If  Y(*)  is  a periodic  autoregression  of  order  (P1»...,Pd)  then 
it  may  be  written  as  (see  Lemma  2) 


(2.3)  L X(t)  + E A'(J)  X(t  - j)  - e'(t) 

J-l 

where  L is  a unit  lower  triangular  matrix,  p satisfies  tie  condition  of 

the  theorem,  and  the  e7(»)  are  uncorrelated  with  a diagonal  covariance 

/ /T  2 2 

matrix,  Ee  (t)  e (t)  ■ D ■ diag  (o, , . . . ,o,)  . Therefore, 

^ ~ 10 

P 

x(t)  + s A(j)  x(t  - j)  - e(t) 
j-l 

where 

A(J)  - L-1A'(j)  J - l,...,p 

e(t)  - L-le'(t)  t-0,±l,.... 


J 
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-1  -T 

and  the  e(*)  are  thus  uncorrelated  and  have  covariance  matrix  E • L DL  , 
which  is  positive  definite  since  D is. 

Conversely,  suppose 

P 

X(t)  + E A(J)  X(t  - j)  - e(t) 
j-1 

where  the  e(*)  are  uncorrelated  with  positive  definite  covariance  matrix 

ryj 

T 

E • Define  the  unique  modified  Cholesky  decomposition  of  E by  E “ LDL  , 

where  L is  unit  lower  triangular  and  D is  positive  definite  diagonal, 

2 2 -1 

D m diag  (0j». . . .rr^)  . Then  L ia  unit  lower  triangular  and  X(*) 

satisfies 

L_1X(t)  + E A'(j)  X(t  - j)  - e'(t) 
j-1 

where 

a'(J)  - L_1A(J)  , j - 1,...,P  , 

and  the  ?/(t)  ■ L * e(t)  , and  are  thus  uncorrelated  with  diagonal  positive 
definite  covariance  matrix,  and  the  theorem  is  proved. 

Note  that  without  too  much  difficulty  we  could  accomodate  a process 
for  which  E is  not  of  full  rank. 

Definition.  A process  Y(*)  is  said  to  be  a covariance  stationary  periodic 
autoregression  if  it  is  a periodic  autoregression  and,  furthermore,  the 
associated  vector  process  X(*)  , (1.1),  is  covariance  stationary. 
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By  a stationary  multiple  autoregression,  we  mean  a process  X(») 
which  obeys  (2.1)  and  admits  to  a solution,  in  the  mean  square  sense,  in 
terms  of  the  "past"  g(#)  , i.e. 


X(t)  • e(t) 


+ S B(k)  e(t  - k) 
k-1 


For  this  to  be  true,  the  zeroes  of  the  determinental  polynomial, 
detfl  + £ A(j)  z^J  must  lie  outside  the  unit  circle.  This  does  not 
translate  into  the  intuitively  expected  conditions  for  the  periodic  auto- 
regression. For  example,  suppose  d * 3 and  p^  «*  p^  - p^  - 1 » then 
for  stationarity  we  require  1^(1)  0^(1)  a^(l)l  < 1 . 

One  of  the  values  of  an  autoregressive  process  is,  of  course,  the 
facility  with  which  one  may  perform  linear  prediction.  It  can  easily  be 
shown,  see  [18],  that  the  least  squares  predictor  of  Y(t  + h)  , for 
positive  integer  h , on  the  basis  of  Y(t),Y(t-l) is 

h-1  pt-rh 

Y(t+h/t)  » - £ a (j)  Y(t  +h  - j/t)  - £ a (j)  Y(t  +h  - j) 
j-1  j=h  t4h 

if  Y(*)  is  a covariance  stationary  periodic  autoregression.  The  one- 
step-ahead  (h  - 1)  prediction  error  is  simply, 

E{Y(t+l/t)  - Y(t+1)}^  - • For  h>  1 one  must  formally 

solve  (see  (3.8)), 

Pt-H, 

Y(t+h)  - - £ a (j)  Y(t  +h  - j)  + e(t  +h) 

J-1  t4h 


1 


“ S fWk>  «<t+h-k) 

k-0  ™ 


am 


j 
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with  g.  (0)  ■ 1 > in  which  case 
t+n 

E{Y(t  +h/t)  - Y(t  +h)]2  - E1  B^OO  aj+h.k  • 

k*0 

A side  benefit  is  that  the  above  method  yields  a solution  to  the 
problem  of  predicting  only  a subset  of  X(t  + l),X(t  +2),...,  in  terms  of 
X(t) ,X(t  - 1),...  , for  an  autoregressive  X(*)  . 
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3.  Parameter  Estimation 

Given  a sample  Y(l) Y(T)  from  a zero  mean  Gaussian  covariance 

stationary  periodic  autoregression  of  order  (p^,...^)  and  period  d , 

2 

we  wish  to  make  inference  about  the  parameters  a^Q)  and  > 

j * l,...,pk  , k = l,...,d  . We  show  that  the  results  of  Mann  and  Wald  [6] 
extend  to  this  case. 

To  simplify  the  notation  we  take  T **  Nd  , where  N is  a natural 
number.  Then  the  problem  is  equivalent  to  having  a sample  X(l) , . . . ,X(N) 
from  a stationary  autoregression,  where  X(*)  and  Y(*)  are  related  by 

(1.1). 

Define 

-1  m 

(3.1)  Vk*v>  **  N S Y(k+dj)  Y(k+v  + dj) 

j«0 

where  m ■ N -[(k  + v)/d]  , for  k = 1,. . . ,d  , v = 0,l,...,T-k-l  . 

Lemma  1.  If  Y(’)  is  a periodically  correlated  Gaussian  process,  then, 
with  R^k.v)  defined  by  (3.1),  as  T -*  » , the  ^(k.v)  converge  almost 
surely  and  in  mean  square  to  R(k,v)  , and 

CD 

(3.2)  N Cov  pMkt.Vj)  , RN(k2,v2)]  -*  S {R(k1,lt2  + du)  R(v1,v2  +du) 

& " u** 

+ R(k1,v2+du)  R(vlt!;2 +du))  . 
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Proof.  By  using  (1.1)  and  (1.2)  and  Gladyshev's  theorem,  we  see  that  this 
is  merely  the  extension  of  Slutsky's  result;  see  [3],  pp.  209,  210. 


In  a periodic  autoregression,  the  tt  and  o are  related  to  the 
covariance  kernel  R in  a modified  Yule-Walker  form; 

Theorem  2.  If  Y(‘)  is  a covariance  stationary  periodic  autoregression 
of  order  (p^,**..?^)  with  covariance  kernel  R(*,*)  , then  for 
k **  1, . . . ,d  , 

pk  2 

(3.3)  R(k,k-v)  + £ C^U)  R(k- j,k-v)  - f)vQak  , v * 0 , 

j=l 

where  we  use  the  Kronecker  delta. 


Proof.  We  must  first  show  that 


rt 

Y(t)  - - £ a (j)  Y(t  - j)  + c(t) 
J-l 


* e(t)  + £ 3t.(k)  e(t-k)  . 
k-1 

This  result  is  iranediately  available  from  Theorem  1.  Therefore,  for  posit? m 
v , Y(t  - v)  is  uncorrelated  with  e(t)  . Multiplying  both  sides  of  (2.2) 
by  Y(t  - v)  and  taking  expected  values  yields  the  theorem. 
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We  are  thus  immediately  led  to  Fisher— consistent  estimators  of  the 
2 

a and  a ; in  (3.3)  replace  the  R by  the  appropriate  R^  and  solve 
the  resulting  linear  equations.  Indeed,  one  can  show  after  some  laborious 
algebra  that  these  estimators  are  the  same  as  those  one  would  obtain  by 
solving  the  appropriately  constrained  multivariate  Yule-Walker  equations. 
Thus,  the  above  suggests  a possible  method  for  finding  the  properties  of 
the  estimators,  but,  since  the  constraints  are  unpleasant,  we  prefer  to 
proceed  directly. 

If  Y(*)  is  a periodically  correlated  process,  denote  by  <9  the 
Fisher  information  matrix,  and  by  J(a,B)  the  element  in  this  matrix 
corresponding  to  the  parameters  a and  £ . 

Lemma  2.  If  Y(*)  is  a Gaussian  covariance  stationary  periodic  auto- 
regression of  order  (p^,...,?^)  , then  the  information  matrix  is  block 
diagonal, 

J(ak(j)’ama))  “ ^ROc-j.k  , 

tKak(J>,^)  **  o 

- 6,„/2<4  • 

for  j * 1| • • • >Pj^  $ i ■ 1, • 


F 
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Proof.  Using  (1.1)  then. 


P 

L X ( t ) S A(1)  X(t  - j)  = r(c) 

1=1 


where 


ho  ’ V1*  • J < k 

Ak.(v)  = O^dv  - k + j)  , v = l,2,...,p 

2 2 

and  p as  in  Theorem  1.  Then,  defining  D = diag  (cr^,  , . . . ,ct, ) , and, 

P 1 i,., 

G(z)  = L £ A(j)  z3  , z = e , 

j-1 

we  have  that  the  spectral  density  matrix  of  the  X(')  process  is 
l'O'O  = G_  1 (z)  DG**(z)/2"  , 


where  the  asterisk  denotes  the  complex  conjugate  transpose.  We  have, 
see  [17],  that 


(3.4) 


4tt 


TT 

P 

af1  f Sf'1  f 

J tr 

-Tf  1 

“ f TT  f. 

(a-)  d 


Let  E,,  denote  the  d-dimensional  square  matrix  with  a one  in  the  (j,k) 

jk 

component  and  zeroes  elsewhere.  Now 


{-^  .. 
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<3,6)  £ " EjkG  * + G*DEkj£  * J<k 

(3.7)  *£l(  . -gVX  • 

k 

Using  the  formula, 

(3.8)  G (z)  - L_1  + S B(k)  zk  , 

k®l 

where  the  B(*)  decay  exponentially  to  zero  because  of  the  assumption  on 
the  zeroes  of  det  ^G(z)j  , we  see  that  the  terms  of  the  form 

tr  £ hi0'1**,0'1'*  1(V'H,) " " ■ 0 

for  all  k,j,m,n  if  v + u ^ 0 , and  for  k>j  or  m>n  if  v + u * 0 
(see  [10]).  Upon  substitution  of  the  appropriate  (3.5),  (3.6)  and,  or,  (3.7) 
into  (3.4),  and  some  straightforward  algebra,  the  lemma  follows. 


If,  with  the  defined  in  (3.1),  we  define  the  a and  a as 

solutions  to  the  normal  equations, 

pk 

(3.9)  R^k.k-v)  + E VJ>  lyk- j,k-v)  - 6^0* 

A A _ 

for  v - 0,...,pk  and  k • l,..,,d  , and  the  vectors  »ak(Pk)]T  » 


then, 
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Theorem  3,  If  Y(l) , . . . ,Y(T)  is  a sample  from  a covariance  stationary 

A 

Gaussian  periodic  autoregression  of  order  (p^ , . . . jp^)  , then  the  a and 

A 

o defined  by  (3.9)  are  almost  surely  consistent  estimators  (T  -*  »)  . And 

«/TT  (ol^  - lc  ■ 1 d have  an  asymptotic  distribution  which  is 

Gaussian  with  mean  zero  and  covariance  matrix  * , where  *9  is  the  appro- 
priate block  of  the  information  matrix  given  in  Lemma  2.  Thus  the  estimators 
are,  in  this  sense,  asymptotically  efficient. 


Proof.  The  consistency  follows  from  Lemma  1 and  Theorem  2.  To  show  the 
asymptotic  distribution,  as  in  [12],  we  have  from  (3.3)  and  (3.9)  for 
k » 1, . . . ,d  and  v = l,...,p  , 

Pk  A Pk 

(3.10)  E Rjk- j,k-v](ok(j)-ak(j)j  - z \U)(Vk-J»k-v)-R(k- J.k-V> 
j°0  j**0 

A 

defining  0^(0)  “ 0^(0)  = 1 . From  Lemma  1 the  random  variables  on  the  right 
of  (3.10)  are  asymptotically  Gaussian  with  mean  zero.  To  find  their  asymptotic 
covariance  matrix,  consider  for  v^  and 


V£  positive. 


14 


Pu  P,. 

1 "2 


N ‘V/V  C°V  •Sl,kl"  vl),Rll(k2‘  j2,k2"  V2*) 

jx=0  j2=0  1 2 


P>.  Pv 
,ll  1'2 


= £ 2 Cs  (j .)  a.  (j,)  (R(k  - j.  ,k  - j,)  R(k  - v.  ,k„  - v?) 

j1=0  j2=0  Li  1 K2  L * 1 1 1 L 1 1 Z 


+ RO^  - ,k. 


- j25 


R(k2  " V2’kl  " ^ 


03 

£ |jl(U1  - j1,k2  - j2 +du)  R(k!  ' Vj^.P-2  “ v2  + du> 


u=l 


+ R(k2  - j2»kj  - J1  + du)  R(k2  - v2 ,k^  - Vj  -idu) 


R(kj_  - ■ 5 2 _!'du>  R^k2  “ v2,kl  " jl  " du) 


+ R(k^  - v1,k2  - J2  - du)  R(k2 


vV  ji 


\6klfk2R<kl-Vl*kl-v2> 
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from  Lemma  1 and  a repeated  application  of  (3.3).  I!ow,  from  Lemma  1 and 
Cramer's  theorem,  ([1]  p.  254)  the  random  variables  on  the  left  of  (3.10) 
have  the  same  joint  asymptotic  distribution  as 


S«**Vk 
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k 

S r[u-  j,k-v](cL  (j)  - cl(J)) 
j=0  u 


v = 1, 

k «*  1, 


*Pk 

,d 


and  the  theorem  is  proved. 


A noteworthy  result  contained  in  the  theorem  is  that,  asymptotically, 

A 

the  , k * l,...,d  , are  independent.  Thus,  when  analyzing  a multivariate 
autoregression,  the  parametrization  in  terms  of  the  , as  opposed  to  the 
A(j)  in  (2.1),  allows  us  to  analyze  each  channel  separately.  In  addition, 
since  the  different  orders  of  the  channels  are  not  necessarily  such  that 

dp  9 Pk  - k + 1 , I:  ■ 1 d , we  can  model  a multivariate  autoregression  of 

order  p (with  A(p)  r 0)  with  fewer  than  d^p  + ^d(d  + l)]/2  parameters, 
and  we  thus  have  a general  methodology  for  systematically  reducing  the  number 
of  parameters  required.  The  number  of  parameters  in  terms  of  the  A(j)  can 
be  sizable;  for  example,  if  d ■ 12  (monthly  data),  to  fit  a zero  order  multi- 
ple autoregression  requires  78  parameters;  a first  order,  222  parameters. 

We  have  thus  overcome  the  difficulty  expressed  by  Whittle  [17],  "...  of 

2 

analyzing  a d-tuple  series  [which]  may  be  said  to  increase  roughly  as  d 
(the  number  of  auto-  and  cross-correlograms  which  must  be  calculated,  and  the 
order  of  the  number  of  parameters  to  be  estimated),  while  the  number  of 
observations  increases  only  as  d ."  It  could  be  said  that  we  have  defined 
multiple  autoregressiono  of  non-integer  orders. 
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4.  Deciding  the  Order 

When  using  autoregressions  to  model  data,  an  important  question  has 

been  what  order  to  use.  To  answer  this  for  periodic  autoregressions,  we 

expand  the  notation  in  (2.1)  and  (2.2)  to  make  explicit  the  orders,  by  using 

2 

Sp  , and  ( Pj. ) . Parzen  [15]  has  introduced  the  CAT  criterion 

to  decide  the  order  of  a fitted  multiple  autoregression, 


CAT  (m) 


1 

tr  [dN  £ 
p=l 


m *=  0,1 


which  measures  the  mean  square  error  (plus  a constant)  of  approximating  the 
infinite  order  autoregressive  transfer  function  by  a fitted  mth  order  auto- 
regression. The  chosen  order  is  that  which  minimizes  an  estimate  of  CAT  (n)  . 

Rewriting  CAT  we  see  that  it  equals 


m d I'  o,2(j,P,,) 

ESS  ~ 

P-1  J=0  Ck(Pk) 


d k 
k-l  j-0 


where  dp  » pfc  - k * 1 , dm  - - k * 1 , k » l,...,d  . 

To  remove  these  constraints  on  the  individual  orders,  and  thus  allow  for 
noninteger  order  multiple  autoregressions,  we  can  define 


d 

(4.1)  PCAT(m)  =>  £ P CAT  (m.  ) 

k-l  lc 


where  m - (m^,...,mj)  , and, 

k f-l”* 

P CAT  (m.  ) - S N £ 

j-o  L p-l  at( p) 


(4.2) 


p) 


1»'V| 


“k  - 0,1,., 


We  see  that  PCAT  (m)  is  minimized  when  each  PCAT  (m^)  is  minimized.  The 

A 

PCAT  (m^)  may  be  estimated  by  replacing  the  a by  the  a in  (3.8)  and  the 
CT^(p)  by  Qj/(N  - p)J  CT2(p)  obtained  in  (3.8).  In  this  form,  we  can  use 
PCAT  to  decide  the  order  of  a periodic  autoregression. 

An  alternative  is  to  use  the  Akaike  criterion  [4) , which  presupposes 
that  the  sample  is  indeed  from  a multiple  autoregression  of  order  m , 

A 2 

AIC  (m)  = log  det  5^  + (2d  m/N)  , 

and  minimize  it  with  respect  to  m . This  function  can  be  rewritten  as 

^ *2  2 

AIC  (m)  » Z log  aAnt  ) + (2d  m/N) 

1=1  lc  K 

with  dm  *»  m^  - k + 1 , k ■ l,...,d  . This  leads  to  the  generalization, 

d 

(4.3)  PA1C  (m)  = Z PAIC  (n^) 

k“l 

where  m * (m^,...,md)  , and, 

(4.4)  PAIC  (m,(<)  = log  0^(0^)  + ^/N  > "V  = °»1 

which  is  a straight  generalization  of  the  univariate  AIC  . This  is  indeed 
approximately  proportional  to, 

— 2 In  (maximum  likelihood  function)  + (2  number  of  parameters)  /II 


agreeing  with  the  definition  given  in  [4] 
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We  note  that,  as  for  scalar  autoregressions,  the  periodically  correlated 
autoregression  allows  us  to  put  a prediction  theoretic  interpretation  on  the 
multivariate  A1C  . If  we  let 

= Nct^C^/CN  - n^)  , k “ 1 d , 

i.e.,  introduce  the  "unbiased"  estimator  of  the  residual  variance,  then 

AIC(m)  - S log  S£(n^)  [l  + n^lf1  + 0(m£  if2)] 
lc=l 

with  dm  ■ m^  - k + 1 , k=l,...,d  . Thus,  asymptotically,  AIC  is  the 

sum  of  the  logs  of  the  estimated  one- step-ahead  prediction  variances  of 
X^(t+1)  given  X(t) ,X(t  - 1) , . . . , plus  that  of  predicting  X2(t+1)  given 
X1(t +1) ,X(t) ,X(t  - 1) , . . . , plus  that  of  predicting  X3(t  +1)  given 
X0(t  +l),X1(t  +l),X(t),X(t  - 1),...,  etc 


I.- 
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5.  Discussion 

Viewing  a multiple  autoregression  as  a periodic  autoregression  clearly 
displays  the  effects  of  prewhitening  each  channel  before  doing  a joint  analysis. 
If  the  prewhitening  is  done  with  an  autoregressive  filter  on  each  channel,  then 
this  can  be  viewed  as,  first,  fitting  a periodic  autoregression  with 
c^(d) ,a^(2d) , . . . , k = l,...,d  , as  the  only  nonzero  coefficients,  and, 
second,  performing  a periodic  autoregression  on  the  residuals.  An  alternative 
approach  would  be  to  include  a subset  regression  option  when  fitting  the 
periodic  autoregression  (such  as  in  [8],  for  example)  and  this  would  obviate 
the  need  for  prewhitening. 

In  analyzing  atmospheric  data,  Jones  and  Brelsford  [5]  achieved  a 
reduction  in  the  number  of  parameters  by  expanding  the  a in  Fourier  series, 

m 

(5.1)  a^Cj)  = s [c  cos  (2rrn  k/d)  + s sin  (2tt  n k/d)]  , 

n=0  J •’n 

taking  m small  (relative  to  d/2  ),  arguing  that  the  a were  slowly  varying 
(with  respect  to  k ) and  periodic  of  period  d . Using  the  asymptotic  dis- 
tribution in  Theorem  3,  one  can  obtain  efficient  estimators  of  the  c's  and 

A 

s's  by  performing  a weighted  regression  of  the  a , in  (3.9),  on  the  c's  and 
s's  (as  in  [111).  Indeed,  these  estimators  are  not  too  different  from  those 
in  [51,  but  this  approach  would  provide  a method  for  systematically  testing 
the  hypothesis  exhibited  in  (5.1). 

The  Gaussian  assumptions  made  in  Section  3 can  clearly  be  relaxed. 

Another  generalization  can  be  achieved  by  considering  periodically  correlated 
q-dimensional  vector  processes  Y(*)  , in  which  case  X(»)  would  be  a 
dq-dimensional  covariance  stationary  time  series  [2l. 
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